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In this paper, we investigate the agegraphic dark energy (ADE) model by includ¬ 
ing the sign-changeable interaction between ADE and dark matter in non-flat 
universe. It is shown that the interaction induces an energy flow of which the 
direction is first from dark matter to ADE and then from ADE to dark matter. 
The phase space analysis is made and the critical points are found, one of which 
is the accelerated scaling attractor solution. So, the coincidence problem can be 
alleviated. Furthermore, we show the evolution of the density parameter D, the 
deceleration parameter q and the equation of state (EoS) parameter Wd of ADE. 
We also find that our model is consistent with the observational data. 
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1. Introduction 


Nowadays it is believed that our universe is accelerating, which is supported by cosmological observations 
[ll-Q . This cosmic acceleration is commonly explained by dark energy with negative pressure. The simplest 
dark energy model is Einstein’s cosmological constant model (ACDM) with EoS wa = —I- Though the 
ACDM can fit the observational data well, the scenario suffers from two serious issues called ’’coincidence 
problem” and ” fine-tuning problem”. Among different candidates for dark energy, the ADE model proposed 
by Cai Q has arisen people’s great interest recently. This model is based on the uncertainty relation of 
quantum mechanics as well as the gravitational effect in general relativity. From quantum fluctuations of 
spacetime, the time parameter t in Minkowski spacetime is not more accurate than 6t = (Karolyhazy 

relation) where 77 is a dimensionless constant of order unity. Base^o_n Karolyhazy relation, the energy 
density of metric fluctuations of the Minkowski spacetime is given by 
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1 
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where tp is the reduced Planck time. On these basis, the energy density of ADE is given by Q 


Pd = 


Sn^rrip 

rp2 
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Here rrip = (SttG) n is a constant parameter and T is chosen to be the age of the universe 


T = 


da 

T lla" 
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where a is the scale factor, H = a/a is the Hubble parameter and a dot denotes the derivative with respect 
to the cosmic time. Furthermore, the interacting ADE has been proposed and investigated Q. It was shown 
that the EoS of interacting ADE can cross the phantom divide Wd = —1 from bottom to top. The interacting 
ADE model also has been extended to the universe with spatial curvature [l^ . 

Since the nature of dark energy and dark matter remains unknown, it will not be possible to derive the 
precise form of the interaction from fundamental theory. One has to discuss it to a phenomenological level. 
The most familiar form of interaction [§ 0 ] between dark energy and dark matter is Q = 3cHp. Here c is 
a coupling constant and positive c means that dark energy decays into dark matter, while negative c means 
dark matter decays into dark energy, and p is taken to be the density of dark matter, dark energy, or the sum 
of them. Obviously, these interactions are either positive or negative and hence can not change their signs. 
However, in Ref.[23j, in a way independent of specific phenomenological form of interaction the authors fitted 




the interaction term Q with observations. It was found that Q was likely to cross the noninteracting line 
(Q = 0), namely the sign of interaction Q changed in the redshift range of 0.45 < 2 ; < 0.9. It is interesting and 
enlightening to consider the possibility that the interaction between dark energy and dark matter changes 
sign during the cosmological evolution. Based on this, Wei [24, 25| proposed a sign-changeable interaction 
as follow 

Q = q{ap + 3l3Hp), (4) 


where a and /3 are both dimensionless constants, q = —aa/a? is the deceleration parameter. It is easy to 
find that the interaction Q can chan ge it s sign when the expansion of our universe changes from deceleration 
{q > 0) to acceleration (g < 0). Wei considered the cosmological evolution of quintessence and phantom 
with this type of interaction and found this type of interaction can bring new features to cosmology. In Refs. 
[26i - l^ , it was also shown that the interaction between dark energy and dark matter can change sign during 
the cosmological evolution. 

On the other hand, it is well known that the flatness of the universe is one of the important predictions 
of conventional inflationary cosmology. The inflation models theoretically produce 12^0 on the order of 10“®. 
But recently Ade et al. [29| (Planck data) found following constraints on 


100II«,o = -4.2t^J(Planck-PWP-fhighL), 

IOOUko = —l.Oi'iJJCPlanck-flensing-l-WP-l-highL). 


These constraints are improved substantially by the addition of BAO data which are 

100n«o = -0.05t^j^(Planck-PWP-hhighL-HBAO), 

IOOUkO = —0.10toJ5(Planck-|-lensing-|-WP-|-highL-|-BAO). 


In addition, since the spatial curvature is degenerate with the parameters of dark energy, it is of great 
importance to study dark energy models with spatial curvature. 

The above discussion motivates us to study ADE model by including the sign-changeable interaction (4) 
between ADE and dark matter in a non-flat universe. The paper is organized as follows: In Sect. 2, we 
make the phase space analysis and discuss the accelerated scaling attractor solution in the interacting ADE 
model. In Sect. 3, we show consistency of the model with the observational data. The conclusions are given 
in Sect. 4. 


2. Dynamical analysis of the interacting ADE model in a non-fiat universe 

We consider the Friedmann-Robertson-Walker (FRW) metric for the non-flat universe as 

ds^ = dt^ — ^-I- (5) 

where k = —1,0,1 is curvature parameter corresponding to open, flat and closed universe, respectively. In 
our scenario, the non-flat FRW universe contains two components, one is ADE component pd and the other 
is dark matter component pm with Wm = 0, ie., the total energy density ptot = Pd + Pm- The corresponding 
Friedmann equation is given by 

We assume that the ADE and dark matter exchange energy through an interaction term Q, namely 


Pd + 3H{1 + Wd)pd = -Q, 

Pm T Pm — Qj 


(7) 

( 8 ) 


where Wd is the EoS parameter of ADE. Since the term ap in (4) is introduced from the dimensional point 
of view [24| . one can remove this term by setting a = 0, and then (4) becomes simply Q = 3l3qHp. In this 
paper, we take the form 

Q = 3PqHpd (9) 

with P being a coupling constant, q is the deceleration parameter 


_ aa ^ 
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( 10 ) 


Obviously, the interaction term can change its sign when the expansion of our universe changes from decel 
eration {q > 0) to acceleration (g < 0). 






In what follows, we investigate the dynamical system by phase space a naly sis and numerically calculate 


the cosmological evolution. For the phase space analysis in ADE, see 
dimensionless variables 

^ Pd I Pm . 

y = J „ o , A = 


e a naiy s 

SIIJ. 
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We introduce the following 


( 11 ) 


Differentiating Eq.(6) with respect to the cosmic time and combing Eqs.(2), (8), (9) and (10), one can easily 
get 

i ?2 1 + 

From Eqs.(ll) and (12), we have the following autonomous system 


x' 

= x(g+l--), 

n 

(13) 

y' 

y , 3 _x^, 

= -5 + (»+2 V‘I. 

(14) 

A' 

= 2gA, 

(15) 

where 

g = -l 

H -l-A+|y2 + ix3 

1 + 1/3x2 

(16) 

and the prime denote the derivative with respect to N = Ino. Erom Eq.(6), we have 



x^ + y^ = 1 + A. 

(17) 


The fractional energy density of the ADE, dark matter and curvature is 

P<^ = 2.2 o = = „2 

Zmim ’ imim ^ ’ 

^ - y 

Ff2a2 

The EoS of ADE could be expressed in terms of the new variables as 

^ , 2x /3(-l - A+|j/2 + i2;3) 

Wd = -^ + t; -;-— - -■ 

3n 1 + |/lx2 


Hd = 


(18) 


(19) 


In the absence of interaction (/3 = 0), from Eq.(19), one can see that Wd is always larger than -1 and cannot 
cross the phantom divide Wd = —1- However, the situation is changed when the interaction term is taken 
into account. In this case (/3 ^ 0), from Eq.(19), it is easy to see that Wd can cross the phantom divide 
provided 


^ ^ /3(-l - A + + Ax^) 

3n 1 + |/3x2 


( 20 ) 


Setting x' = j/' = A' = 0 in Eqs. (13), (14) and (15), we can obtain the physically meaningful critical 
points (xc, j/c) Ac) of the autonomous system. The critical points and their properties are given in the Table 
1. The first critical point, (0,1,0) is unstable and corresponds to a matter dominated phase. The other 
critical point, (x*,y*,0) is the attractor, x* satisfies 


— —^x^ + (3 — 3/3)x^ + — —3 = 0 
n n 


( 21 ) 


and y* = y^l — x2. The deceleration parameter and the EoS of ADE in the attractor are equal respectively 
< 7 * = —1 + ^ and Wd* = —1 + /I + (| — /3)^. Eor example, consider the case when n = 3.5, /3 = —0.3. 
We can find that in this model y* « —0.75 and Wd* ss —1.05. In the absence of interaction (/3 = 0), from 
Eq.(21), one can get x* = 1, which corresponds to ADE dominated phase. Eq.(21) can be also expressed as 


r(3/3-2)x^ + 2l(--l) = l-x2. (22) 

n 

Since xj < 1, from Eq.(22) we can find /3 < 0 for n > 1. It means that the interaction term (9) can change 
its sign from Q < 0 to Q > 0 when the expansion of our universe changes from deceleration {q > 0) to 
acceleration (g < 0). If the attractor (x*, y*, 0) is the scaling solution, the coincidence problem is alleviated. 

















Table 1: Location of the critical points of the autonomous system of Eqs. (13-15), their stability and 
dynamical behavior of the Universe at those points. 


(^C5 Vc-j ^c) 

Stability 

q 

Wd 

coordinates 

character 



(0,1,0) 

unstable 

1/2 

-1-^/2 

(x,,y*,0) 

attractor 

g* < 0 

Wd* 



Figure 1: The phase plane for n = 3.5 and P = —0.3. The red dot stands for the late-time attractor with 
cc* ~ 0.89, j/* ~ 0.46 (Ud* ss 0.79, ftm* ~ 0.21). 


Scaling solutions are characterized by a constant dark matter to dark energy ration r = jx^. We can find 
that 

/ = 3 / 3 ( 7 -(3-— )r. (23) 

n 

The scaling solutions mean r' = 0, which results in r* = — n)/(3n — 2a::*). So the attractor (a;*, y*, 0) is 

just the scaling solution for (/3 < 0). Thus the coincidence problem can be alleviated because, regardless of 
the initial conditions, the system evolves toward a final state where the ratio of dark matter to dark energy 
stays constant. In Fig.l, it is easy to find that r* 0.27. In addition, since x* < 1, we can see ( 7 * < 0 for 
n > 1. So the critical point (x*,y*,0) is the accelerated scaling attractor solutions. 

In Fig.2, we show the evolution of and Vim with respect to N = In a. We fix n = 0.35 and let /3 vary. 
From Fig.2, one can see that increasing the value of |/3| results in lager Vtm and smaller Vld in the late time. 

Fig.3 shows the evolution of the deceleration parameter q with respect to N = In o for different /3 and the 
fixed n. From Fig.3, we can see that the universe has a transition from deceleration {q > 0) to acceleration 
{q < 0). In addition, the model parameter may affect the time of onset of the acceleration. For the larger 
\(5\ the acceleration sets in earlier. 

In Fig.4, the evolution of FoS of ADF Wd is shown for different /3 and the fixed n. We find that Wd is 
more likely to cross the phantom divide Wd = —V from top to bottom with the increasing of the |/3| in our 
model. Whereas in Ref.[9|, the FoS of interacting ADF can cross the phantom divide from bottom to top. 

3. Observational data 

Since there is mo fundamental theory which can be used to select a specific interacting dark energy 
model, any interacting dark energy model will necessarily be phenomenological. There are two criterions to 
determine whether the model is correct and feasible. One is to examine whether the interacting model can 
lead to the accelerated scaling attractor solutions, which alleviate the coincidence problem, the other is the 
observations. In Sect. 2, it has been shown that the accelerated scaling attractor solutions do exist in our 
model. In the present section, we will consider the Union2.1 Type la supernova dataset (SNela) [s^ , which 
contains 580 data points. We will check the interacting ADE model via the observational data. 





















Figure 2: The evolution of fid (solid line) and ftm (dot line) for n = 3.5 and /3 as indicated. 
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Figure 3: The evolution of deceleration parameter q for n = 3.5 and /3 as indicated. 
































Figure 4: The evolution of the EoS of ADE Wd for n = 3.5 and f3 as indicated. 
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Figure 5: The observed 580 SNela distance modulus along with the theoretically predicted curves (red solid 
line) in the model with n = 3.5 and /3 = —0.3, where we take a priori that current dimensionless Flubble 
parameter h = 0.70. 


The 580 data points of SNela are given in terms of the distance modulus fiobs{zi)- On the other hand, 
the theoretical distance modulus is defined as 

fithizi) = 5\ogiQD L{zi) + Ho, (24) 

where the zero offset fio depends on Hq (or h) as 

Ho = 5 logio + 25 = -5 logio h + 42.38 (25) 

and h is the Hubble constant Hq in units of 100 km s“^Mpc“^. The theoretically predicted value Dl{z) in 
the context of a given model H{z;p) can be described by 

where p denotes the model parameters, = ^/{HoaQ) and Sinn(a;) = sin(a:),x,sinh(a::) for respectively a 
spatially closed > 0, flat = 0 and open < 0 universe. 

Different models result in different theoretical distance modulus Hth- We can judge the plausibility of 
an cosmological model by comparing Hth and the observational value of Hobs- In order to see whether 
the theoretical model corresponds to the observational data, we must have the value of H, which can be 
obtained through the autonomous system (13-15). From Fig.5, one can see that our model is consistent with 
the observational data. 

4. Conclusions 

In the present paper, we investigate ADE model by including the sign-changeable interaction {Q = 
ipqHpd) between ADE and dark matter in non-flat universe. Using the phase-plane analysis, the dynamical 





















behavior of the model has been studied. It was found that the accelerated scaling attractor solution did 
exist in the model, which can help to alleviate the coincidence problem. Moreover, we find that the the 
condition /3 < 0 is necessary for the existence of the scaling solution. It means that the interaction Q can 
change its sign from Q < 0 to Q > 0 when the expansion of our universe changes from deceleration (q > 0) 
to acceleration (g < 0), which is opposite to the case in Ref. . This indicates that at first dark matter 
decays to ADE, and then ADE decays to dark matter. We also show numerically that the coupling constant 
/3 plays an important role in the evolution of the universe. It is easy to see from Figs.2-4 that for the fixed 
n, increasing the value of |/3| can bring the features as follows: larger 17^ is obtained in the late time, the 
acceleration sets in earlier and Wd is more likely to cross the phantom divide from top to bottom. 

Next, using 580 SNela data, we checked the interacting ADE model. The model has given a series of 
reasonable pictures of the cosmic evolution and it is consistent with the observational data. Our work implies 
that we should pay more attention to a sign-changeable interaction between dark sectors. 
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